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Comment on "Dynamic behavior of anisotropic 
non-equilibrium driving lattice gases" 

In a recent Letter Albano and Saracco (AS) study 
the dynamic behavior of some anisotropic driven lattice 
gases by Monte Carlo (MC) simulations [2j. They check 
the scaling of the order parameter OP defined in their 
Eq. (4) and, by data collapse, determine the relevant 
scaling exponents. Then, they formulate scaling Ansatze 
(see Eqs. (5) and (8) of AS) in order to relate the mea- 
sured scaling exponents with the critical exponents an- 
alytically computed within different field-theoretical ap- 
proaches (see Refs. 0, 0)- Comparing MC results and 
theoretical predictions they conclude that the universal- 
ity class of the phase transition in the driven lattice gases 
considered in the paper is well-described by the field the- 
ory proposed in Ref. , ruling out that of Ref . |s| • In this 
Comment we point out that the Ansatze (5) and (8) do 
not take properly into account the strongly anisotropic 
nature of the phase transition, by implicitly assuming 
z = zj_ = Z||. As a consequence, at variance with AS's 
claims, MC data are not conclusive on the issue of the 
universality class. 

Let us outline the argument. During the relaxation of a 
spin system from an initial uncorrelated state with small 
magnetization mo, the /c-th moment of the magnetization 
is expected to scale according to [2j M^\t, t,L, mo) = 
b~ kl3 / v f k {b' z t,b 1 / w ,b 1 / v T,b- 1 L,b Xo m ), where x 
is related to the initial-slip exponent. The gen- 
eralization of this expression to a system that is 
strongly anisotropic and has a conserved order pa- 
rameter requires some care. Using the notations 
of Ref. H, one expects M^ k \t, r, Ln,L±,mo) = 
b~ k P/ v w f k (b- z « t, b x l v w r, b- x L\\ , b- v ^l v W L ± ,b xa m ), 
where M can be naturally identified with some suitable 
order parameter, for instance with the excess density 
as in AS. Since mo is approximately zero for the initial 
configurations, taking into account that the OP defined 
by AS is expected to scale as V Af( 2 ) lj , one obtains Eq. 
(5) of AS with the identification z = Zn . Correctly, in 
Table I, AS report the predictions z = 4/3 and z ss 1.998 
for the two models they consider (see Eqs. (1) and (2) 
of AS), which refer indeed to z\\. Formulae (6) and (7) 
follow from Eq. (5) and therefore they are correct only 
if we identify z = z\\. In Eq. (8), AS generalize the 
scaling behavior of OP to the case of the relaxation 
from the ground state configuration. Without attemping 
a derivation of their equation, we simply note that, 
as it stands, it is not compatible with Eq. (5) of AS: 
The exponent z appearing in Eq. (8) must be identified 
with z± — v\\Z\\/vx and not with z». Indeed, whatever 
the initial condition is, we expect that after a time 
interval long enough to allow the relaxation of the initial 
condition, OP will display a characteristic long-time 
behavior, that should be captured by any reasonable 


scaling form. Now, from Eq. (5) (rewritten as before) 
it is easy to derive that, in the thermodynamic limit, 
OP(t,r) = t-^/ u \\ z \\f FT >(t 1/v ^r), while from Eq. (8) 
OP(t,r) = t-P/ u ^ z f GS (t 1 / u ^ z T). These two equations 
give the same long-time limit only if in the second 
case we identify z with zj_, i.e. OP(t,T, L», L±) = 
b-P/^OP{b- z ^t, b^r, b-^i/^L^b^Lx). Equations 
(9) and (10) follow from Eq. (8) and therefore here z 
must be interpreted as z±. As a consequence, some of 
the theoretical predictions for the quantities measured by 
AS and reported in Tab. I of Ref. pj, change according 
to the following table 


Model 


cx(+) 

Eq. (1) in Ref. [1] 

1/4 

1/4 

Eq. (2) in Ref. [lj 

» 0.13 

re 0.27 


Comparing these corrected predictions with the numer- 
ical results reported in Tab. I of AS, we conclude that 
the value of {$jv\_z is in agreement with Eq. (1), which is 
exactly the opposite conclusion with respect to AS, while 
the numerical estimate of c± differs from the predictions 
of both field-theoretical models considered in AS's paper. 
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